Anomalous scaling of a passive scalar advected by the Navier— Stokes velocity field: 

Two-loop approximation 
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The field theoretic renormalization group and operator product expansion are applied to the 
model of a passive scalar quantity advected by a non-Gaussian velocity field with finite correlation 
time. The velocity is governed by the Navier-Stokes equation, subject to an external random 
stirring force with the correlation function oc S{t — t')k'^~'^~^'^ . It is shown that the scalar field is 
. intermittent already for small e, its structure functions display anomalous scaling behavior, and the 

' corresponding exponents can be systematically calculated as series in e. The practical calculation 

, is accomplished to order (two- loop approximation), including anisotropic sectors. Like for the 

' well-known Kraichnan's rapid-change model, the anomalous scaling results from the existence in the 

model of composite fields (operators) with negative scaling dimensions, identified with the anomalous 
r-{ , exponents. Thus the mechanism of the origin of anomalous scaling appears similar for the Gaussian 

.^fH ' model with zero correlation time and non-Gaussian model with finite correlation time. It should be 

I emphasized that, in contrast to Gaussian velocity ensembles with finite correlation time, the model 

■ and the perturbation theory discussed here are manifestly Galilean covariant. The relevance of these 

results for the real passive advection, comparison with the Gaussian models and experiments are 
briefiy discussed. 



PACS number(s); 47.27.-i, 47.10. +g, 05.10.Cc 



I. INTRODUCTION 



In recent years, considerable progress has been achieved in the understanding of intermittency and anomalous 
^ • scaling of fluid turbulence. Both the natural and numerical experiments suggest that the deviation from the clas- 
\ sical Kolmogorov theory [1-3] is even more strongly pronounced for a passively advected scalar field 9{x) = 9(t,yi) 
> (temperature, entropy, density of an impurity etc) than for the velocity field itself; see e.g. Ref. [3-6] and literature 
' cited therein. At the same time, the problem of passive advection appears more easily tractable theoretically: even 
OO simplified models describing the advection by a "synthetic" velocity field v(a;) = {vi{x)} with a given Gaussian statis- 
tics reproduce many of the anomalous features of genuine turbulent heat or mass transport observed in experiments. 
Therefore, the problem of passive scalar advection, being of practical importance in itself, may also be viewed as a 
starting point in studying intermittency and anomalous scaling in the turbulence on the whole. 

The issue of interest is, in particular, the behavior of the equal-time structure functions of the scalar field 



Ch' 5„(r) = ([0(^,x)-^^(^,x')]"), r = x-x', r = |r|. (1.1) 

> ■ 

The concept of anomalous scaling implies a power-law behavior of the functions (1.1) in the inertial-convective range 
of scales, Sn oc r^", with a nonlinear dependence of the exponents ^„ on n; see e.g. Refs. [1-6]. 

The crucial breakthrough in theoretical research is related to a simple model of a passive scalar quantity advected 
by a random Gaussian field, white in time and self-similar in space, known as the Kraichnan's rapid-change model [7]. 
There, for the first time the existence of anomalous scaling was established on the basis of a microscopic model [8] 
and the corresponding anomalous exponents were calculated within controlled approximations [9-12] and a systematic 
perturbation expansion in a formal small parameter [13]. Detailed review of the recent theoretical research on the 
passive scalar problem and the bibliography can be found in Ref. [14]. 

In the "zero- mode approach," developed in [9-12] (see also [14]), nontrivial anomalous exponents are related to 
the zero modes (unforced solutions) of the closed exact differential equations satisfied by the equal-time correlation 
functions. From the field theoretic viewpoint, this is a realization of the well-known idea of self-consistent (bootstrap) 
equations, which involve skeleton diagrams with dressed lines and dropped bare terms. Owing to very special features 
of the rapid-change models (linearity in the passive field and time decorrelation of the advecting field) such equations 
are exactly given by one-loop approximations, and the resulting equations in the coordinate space are differential 
(and not integral or integro-differential as in the case of a general field theory). In this sense, the model is "exactly 
solvable." Furthermore, in contrast to the case of nonzero correlation time, closed equations are obtained for the 
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equal-time correlations, which arc Galilean invariant and, therefore, not affected by the so-called "sweeping effects" 
that would obscure the relevant physical interactions. 

The second systematic analytical approach to the rapid-change model, proposed in papers [13], is based on the field 
theoretic renormalization group (RG) and operator product expansion (OPE). There, anomalous scaling emerges as a 
consequence of the existence in the model of composite fields with negative scaling dimensions ( "dangerous composite 
operators"), identified with the anomalous exponents. This allows one to give alternative derivation of the anomalous 
scaling, to construct a systematic perturbation expansion for the anomalous exponents, analogous to the famous e 
expansion in the RG theory of critical behavior (see the monographs [15,16] and references therein), and to calculate 
the exponents to the second [13,17] and third [18,19] orders. 

The two approaches complement each other well: the zero-mode technique allows for exact (nonperturbative) 
solutions for the anomalous exponents related to second-order correlation functions [10,20,21] (they are nontrivial 
for passive vector fields or anisotropic sectors for scalar fields), while the RG approach form the basis for systematic 
perturbative calculations of the higher-order anomalous exponents. For anisotropic velocity ensembles or/and passively 
advected vector fields, as well as for passive advection of extended objects (polymers or membranes), where the 
calculations become rather involved, all the existing results for higher-order correlation functions were derived only 
by means of the RG approach and only to the leading order in e [22-25] . 

From a more physical point of view, zero modes can be interpreted as statistical conservation laws in the dynamics 
of particle clusters [26]. The concept of statistical conservation laws appears rather general, being also confirmed by 
numerical simulations of Refs. [27,28] , where the passive advection in the two-dimensional Navier-Stokes (NS) velocity 
field [27] and a shell model of a passive scalar [28] were studied. This observation is rather intriguing because in those 
models no closed equations for equal-time quantities can be derived due to the fact that the advecting velocity has a 
finite correlation time (for a passive field advected by a velocity with given statistics, closed equations can be derived 
only for different-time correlation functions, and they involve infinite diagrammatic series). 

One may thus conclude that breaking the artificial assumptions of the time decorrelation and Gaussianity of the 
velocity field is the crucial point. 

Besides the calculational efficiency, an important advantage of the RG approach is its relative universality: it is 
not bound to the aforementioned "solvability" of the rapid-change model and can also be applied to the case of finite 
correlation time or non-Gaussian advecting field. In Refs. [29,30] (see also [31] for the case of compressible fiow and 
[32] for a passive vector field) the RG and OPE were applied to the problem of a passive scalar advected by a Gaussian 
self-similar velocity with finite (and not small) correlation time. The energy spectrum of the velocity in the inertial 
range has the form £{k) oc k^^"^^ , while the correlation time at the momentum k scales as t{k) oc fc^^+''. It was shown 
that, depending on the values of the exponents e and ry, the model reveals various types of inertial-range scaling 
regimes with nontrivial anomalous exponents, which were explicitly derived to the first [29,31] and second [30] orders 
of the double expansion in e and rj. The most interesting case is ry = e, when the exponents can be nonuniversal 
through the dependence on the ratio of the velocity correlation time and the turnover time of the passive scalar. 

Earlier, a similar model was proposed and studied in detail (using numerical simulations, in two dimensions) in 
[33] . Various aspects of the transport and dispersion of particles in random Gaussian self-similar velocity fields with 
finite correlation time were also studied in Refs. [34-41]. 

As was pointed out in Ref. [33], the Gaussian model with finite correlation time suffers from the lack of Galilean 
invariance and therefore misrepresents the self- advection of turbulent eddies. It is well known that the different-time 
correlations of the Eulerian velocity field are not self-similar, as a result of these "sweeping effects," and depend 
substantially on the integral scale. It would be much more appropriate to impose the scaling relations for £{k) and 
t{k) in the Lagrangian frame, but this is embarrassing due to the daunting task of relating Eulerian and Lagrangian 
statistics for a flow with a finite correlation time.^ In the RG and OPE formalism, the sweeping by the large-scale 
eddies is related to the contributions of the composite operators built solely of the velocity field v(a;) and its temporal 
derivatives, as discussed in detail in Refs. [16,43-46] for the case of the stochastic NS equation. In the Gaussian model 
those operators become dangerous (that is, their scaling dimensions become negative) for e > 1/2, which gives rise 
to strong infrared divergences in the correlation functions [29]. This means that the sweeping effects, negligible for 
small e's, become important for e > 1/2. In a Galilean-invariant model, such operators give no contribution to the 
quantities like (1.1), as explained in [43 46] for the NS case. In the Gaussian case, these IR divergences persist in the 
structure functions, which provides not only an upper bound for the reliability of the e-^ expansion, but also a natural 



^In this connection, it should be noted that, due to the time decorrelation, in the rapid-change model there is no problem in 

relating Eulerian and Lagrangian statistics of the velocity field: they are identical. This allows one to perform very accurate 
numerical simulations in the Lagrangian frame; see Refs. [42]. 
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bound for the validity of the Gaussian model itself (which exchides, in particular, the most realistic Kolmogorov's 
value e = 4/3 and its vicinity). These conclusions agree with the nonperturbativc analysis of Ref. [41], where the 
value of e = 1/2 was reported as the threshold between two qualitatively different regimes for a Lagrangian particle 
advected by a Gaussian velocity ensemble. The same threshold value of e = 1/2 was obtained earlier in Refs. [37] for 
a two-dimensional strongly anisotropic model. 

In this paper, we shall study the case of a passive scalar field, advected by a non-Gaussian velocity field governed 
by the stochastic NS equation. To be precise, the advection-diffusion equation for the scalar field has the form 

Vte = Kod^e + f, Vt = dt + ividi), (1.2) 

where dt = d/dt, di = d/dxi, Vt is the Lagrangian derivative, kq is the thermal conductivity or molecular diffusivity, 
d"^ is the Laplace operator, and / = f{x) is an artificial Gaussian random noise with zero mean and correlation 
function 

{f{x)f{x'))^6{t-t')C{r), r = x-x'. (1.3) 

The detailed form of the function C is unessential; it is only important that C decreases rapidly for r ^ L, where L 
is some integral scale. The noise maintains the steady state of the system and, if C depends on the vector r and not 
only on its modulus r = |r|, is a source of large-scale anisotropy. With no loss of generality it can be assumed that 
the function C is dimensionless. The absence of the time correlation in (1.3) is also unessential; in a more realistic 
formulation, the noise can be replaced by an imposed constant gradient of the scalar field; see Refs. [11,12,29,31,33]. 

The transverse (divergence-free, due to the incompressibility condition = 0) velocity field satisfies the NS 
equation with a random driving force 

VtVi = uod^Vi-diV + fi, (1.4) 

where P and fi are the pressure and the transverse random force per unit mass (all these quantities depend on x). 
We assume for / a Gaussian distribution with zero mean and correlation function 

(/i(x)/,-(x')) = / rfkP,,.(k)rf/(fc)exp[ik(x-x')], (1.5) 

[^T^) Jk>m 

where -Pij(k) = Sij — kikj/k"^ is the transverse projector, df{k) is some function of A; = |k] and model parameters, and 
d is the dimension of the x space. The momentum m = l/£, the reciprocal of another integral scale i, provides IR 
regularization. Its precise form is unessential, the sharp cutoff is the most convenient choice from purely calculational 
reasons. For definiteness, in what follows it is always assumed that L:^ £, that is, the largest scale in the problem is 
the integral scale related to the scalar noise; it will be set to infinity whenever possible. 

The standard RG formalism is applicable to the problem (1.4), (1.5) if the correlation function of the random force 
is chosen in the power form 

df{k) = Dok^-''-^^ (1.6) 

where Dq > is the positive amplitude factor and the exponent < £ < 2 plays the role analogous to that played by 
4 — d in the RG theory of critical behavior [15,16]. The form (1.6) is widely used in the RG theory of turbulence since 
the pioneering works [47-50]. The most realistic value of the exponent is e = 2: with an appropriate choice of the 
amplitude, the function (1.6) for e ^ 2 turns to the delta function, df{k) cx (5(k), which corresponds to the injection 
of energy to the system owing to interaction with the largest turbulent eddies; for a more detailed discussion see Refs. 
[16,44-46]. 

The results of the RG analysis of the model (1.4)-(1.6) are reliable and internally consistent for small e, while the 
possibility of their extrapolation to the real value e = 2 and thus their relevance for the real fluid turbulence is far 
from obvious; see e.g. Ref. [46] for a recent discussion. We shall not discuss this important problem in detail and 
restrict ourselves with a few remarks which will be relevant in what follows. 

The time decorrelation of the random force guarantees that the full stochastic problem (1.2)-(1.5) is Galilean 
invariant for all values of the model paranic;ters, including Dq and e. As a consequence, the ordinary perturbation 
theory for the model (that is, the expansion in the nonlinearities, or, equivalently, in Dq from (1.6)), is manifestly 
Galilean covariant: all the exact relations between the correlation functions imposed by the Galilean symmetry (Ward 
identities) are satisfied order by order. The rcnormalization procedure does not violate the Galilc^an symmetry, so 
that the improved perturbation expansion, obtained with the aid of RG and OPE, remains covariant. This means, 
in particular, that the Galilean invariant quantities, for example, the equal-time structure functions (1.1), are not 
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affected by the sweeping (here, the latter becomes important for £ > 3/2 [43]). More formally, the contributions of the 
"dangerous" operators built of the velocity field and its temporal derivatives, do not appear in the OPE for invariant 
correlation functions; see Refs. [16,43-46] for detailed discussion. 

This means that the scaling relations, obtained for Galilean invariant quantities for small e, in model (1.2)-(1.6) 
can be extrapolated beyond the threshold £ = 3/2 despite the fact that the sweeping becomes important there. More 
physically, this means that, in contrast to the Gaussian model, the relative motion of the fluid or impurity particles 
in the inertial-convective range of scales is not affected by overall sweeping by the large-scale eddies. Indeed, the 
most recent numerical simulations of the model (1.4)-(1.6) have shown that the scaling relations, obtained by the RG 
analysis for the structure functions, remain valid for e as high as £ = 7/4 [51] (sec also an earlier work [52]). 

For small £, critical dimensions of all composite operators in the model (1.4)-(1.6) are positive. As a result, the 
scaling behavior of the velocity correlation functions is not anomalous, in the sense that they have a finite limit for 
a = l/m — > oo, and the corresponding scaling exponents are multiples of a single quantity (critical dimension of the 
velocity field). 

However, the same numerical simulations [51,52] suggest that, as e increases, the behavior of the model (1.4)- 
(1.6) undergoes a qualitative changeover and the scaling becomes anomalous, in the sense that the exponents of the 
structure functions become different from the results of naive extrapolation of the small-£ prediction and, probably, 
become independent of e. Wo shall return to this important issue in the Conclusion, and in the bulk of the paper we 
shall concentrate on the behavior of the passive scalar field in the model (1.2)-(1.6), which appears highly nontrivial. 

We will show that, already for infinitesimal values of e, when the velocity statistics is not yet intermittent, the 
scalar field, advected by such velocity ensemble, displays anomalous scaling behavior. The corresponding anomalous 
exponents can be calculated within a systematic perturbation expansion, as series in e. 

The plan of the paper is the following. 

Detailed description of the model has already been given above. In Sec. II wc give the field theoretic formulation 
of the original stochastic problem and present the corresponding diagrammatic technique. In Sec. Ill wc analyze UV 
divergences in the model, establish its multiplicative renormalizability, and derive the corresponding RG equations. 
In Sec. IV we show that the RG equations of our model have the only IR attractive fixed point in the physical range 
of parameters; its coordinates are calculated to the second order of the e expansion. Existence of such fixed point 
means that the correlation functions of our model in the IR range (l/r~fc~TO<cA) exhibit scaling behavior with 
certain critical dimensions Ap oi all the fields and parameters F. This determines the dependence of the correlation 
functions on the argument Ar (but not yet on mr). In general, the dimensions are calculated as series in £, but for 
some basic quantities (velocity and scalar fields, their powers, and frequency) they are found exactly. 

The key role in the following is played by the composite operators built of the gradients of the scalar field. They are 
introduced in Sec. V, and the corresponding dimensions Ai? are given to the first order in e (one-loop approximation). 
In Sec. VI we introduce the operator-product expansion and demonstrate its relevance to the issue of inertial-range 
anomalous scaling. We show that the critical dimensions of aforementioned composite operators can be identified with 
the anomalous exponents, which describe the dependence of the correlation functions on the argument mr. The scalar 
operators are "dangerous" (Air < 0), which implies anomalous scaling (singular dependence on mr and divergence for 
mr 0). The anomalous exponents of anisotropic contributions are determined by the critical dimensions of tensor 
composite fields and thus they also exhibit nontrivial scaling behavior. 

The largest Sec. VII is devoted to the calculation of the anomalous exponents (critical dimensions of the composite 
operators built of the gradients of the scalar field) to the order (two- loop approximation). The results look rather 
cumbersome and are presented in a separate section. Sec. VIII. The discussion of the results, their relevance for the 
real passive advection, comparison with the Gaussian models and experiments are given in Sec. IX. 



The field theoretic formulation and renormalization of the problem (1.2)-(1.6) is discussed in detail in Ref. [16,44,45]; 
below we confine ourselves to only the necessary information. 

According to the general theorem [53], stochastic problem (1.2)-(1.6) is equivalent to the field theoretic model of 
the doubled set of fields $ = {v', v, 9' , 6} with action functional 



II. FIELD THEORETIC FORMULATION 



v', v) + 9'De9'/2 + 9' [-Vt + Kod^] 9, 



(2.1) 



where 



5„(v', v) = + v' [-V( + vod""] 



V 



(2.2) 
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is the action functional for the stochastic problem (1.4)-(1.6), Dg and Dy are the correlation functions (1.3) and (1.5) 
of the random forces / and /j, respectively, and all the required integrations over x = {t, x} and summations over the 
vector indices are understood, for example, 

v'{vd)v = dt dx v[{vjdj)v[. 



The auxiliary vector field is also transverse, div[ = 0, which allows to omit the pressure term on the right-hand side 
of Eq. (2.2), as becomes evident after the integration by parts: 



J dtj dx v'idiV = - j df j dx V{div[) = 0. 



Of course, this does not mean that the pressure contribution can simply be neglected: the field v' acts as the transverse 
projector and selects the transverse part of the expressions to which it is contracted in Eq. (2.2). 

Formulation (2.1), (2.2) means that statistical averages of random quantities in the original stochastic problem can 
be represented as functional averages with the weight expS'(<l?), and the generating functional of total [G(^)] and 
connected [W(A)] correlation functions of the problem are represented by the functional integral 

G{A) = exp W{A) = j exp[S'($) + (2.3) 

with arbitrary sources A = {A'"' , A'" , A^' , A^} in the linear form A^ = / dxA^{x)^{x). 

The model (2.2) corresponds to a standard Fcynman diagrammatic technique; the bare propagators (lines in the 
diagrams) in the frequency momentum (lu k) representation have the forms 

(viv'j)^ = {v',vj)l = {-ioj + uok^)-' PM, (viVj)^ = (a;2 + u^.k^' d}{k)P,^{\,), {v[v'^)^ = (2.4) 

with df{k) from Eq. (1.6). The interaction in (2.2) corresponds to the triple vertex —v'{vd)v = vlVijsVjVs/2 with 
vertex factor 

Vijs = i(kj5is + ks5ij), (2.5) 
where k is the momentum argument of the field v' . The full problem (2.1) involves additional propagators 

{ee')^ = {-\u; + Koe)-\ {e6)^ = c{\^){uj^ + i4k'')-\ {o'e')^ = Q, (2.6) 

where C(k) is the Fourier transform of the function C(r) from (1.3); the additional vertex factor in 9'{vd)6 = O'ViViO 
has the form 

Vi = iki, (2.7) 

where k is the momentum of the field 6'. 



III. RENORMALIZATION AND RG EQUATIONS 

The analysis of UV divergences is based on the analysis of canonical dimensions. In contrast to static models, 
dynamical models of the type (2.1), (2.2) have two scales, i.e., the canonical dimension of some quantity F (a field or 
a parameter in the action functional) is described by two numbers, the momentum dimension dp and the frequency 

dimension d'^. They are determined so that [F] ~ [L]^'^''[T]^'^p , where L is the length scale and T is the time scale. 
The dimensions are found from the obvious normalization conditions = —d^ = 1, = = 0, d^ = d^ = 0, 
d!^ = —d,f = 1, and from the requirement that each term of the action functional be dimcnsionlcss (with respect to the 
momentum and frequency dimensions separately). Then, based on dp and dp, one can introduce the total canonical 
dimension dp = dp + 2dp (in the free theory, dt oc 9^), which plays in the theory of renormalization of dynamical 
models the same role as the conventional (momentum) dimension does in static problems. 

The canonical dimensions for the problem (2.1), (2.2) are summarized in Table I, where we introduced the new 
parameters ( "coupling constants" or "charges" ) 

9o = Do/vl, uo = Ko/vo (3.1) 
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instead of Dq and kq and included the dimensions of renormalized parameters, which will appear later on. The 
dimensionless ratio uq has the meaning of the reciprocal of the Prandtl number. From Table I it follows that the 
model (2.1), (2.2) is logarithmic (the coupling constant go is dimensionless) at e = 0, and the UV divergences have 
the form of the poles in e in the correlation functions of the fields 

The total canonical dimension of an arbitrary 1-irreducible correlation function F = ($ • • • $)i-ir is given by the 
relation dr = d'^ + 2(ip = {d+2) — N^d^, where = {N^, N^, Ng',Ng} are the numbers of corresponding fields 
entering into the function F, and the summation over all types of the fields is implied. The total dimension dr plays 
the role of the formal index of the UV divergence: superficial UV divergences, whose removal requires counterterms, 
can be present only in those functions F for which dr is a non-negative integer. Analysis of the divergences in our 
model should be augmented by the following considerations: 

(i) For any model with the Martin-Siggia-Rose-type action, that is, the action of the form (2.1), (2.2), all the 
1-irrcduciblc hmctions with TVv' = Ng' = contain closed circuits of retarded propagators and vanish. 

(ii) If for some reason a number of external momenta occurs as an overall factor in all the diagrams of a given 
Green function, the real index of divergence rfp is smaller than dr by the corresponding number of unities. The 
correlation function requires counterterms only if dp is a non-negative integer. In our model, the derivative d at the 
vertices v'{vd)v and 0'{vd)0 can be moved onto the fields v' and 0' using the integration by parts, by virtue of the 
transversality of the field v. This decreases the real index of divergence: d'^ = dr — N^' — Ng' — Ng, and the fields v', 
0' and enter the counterterms only in the form of the derivatives, dv' and so on. 

(iii) From the explicit form of the vertex and bare propagators it follows that A^^i' — Ng = 2A^o for any 1-irreducible 
function, where TVq > is the total number of bare propagators {00) a entering into the hmction (obviously, no 
function with A^o < can be constructed). Therefore, the difference Ngi — Ng is an even non-negative integer for any 
nonvanishing function. This is a consequence of the linearity of the original stochastic equation (1.2) in the field 0. 

(iv) Galilean symmetry of our problem requires that the counterterms to the action be invariant. In particular, the 
monomials v'dtv, v'{vd)v, 0'dt0, and 0'{vd)0 can appear only in the form of covariant derivatives u'Vf?; and 0'\7t0- 

From Table I we find dr = (d + 2) - (d - l)Ny' - Ny+Ng-{d+ l)Ng> and dj, = (d + 2) - dN,,> - Ny - (d + 2)Ng'. 
Bearing in mind that Ng> > Ng we find that superficial divergences can only be present in the 1-irreducible functions 
{v'v)i-ir and {0'0)i-ir, and the corresponding counterterms reduce to the forms v'd^v and 0'd'^0. The monomials 
v'dfV and 0'dt0 do not contain spatial derivatives and therefore they are ruled out by the property (ii). Then the 
property (iv) excludes the monomials v'{vd)v and 0'{vd)0 (allowed by dimensional considerations). 

In the special case d = 2 a new UV divergence appears in the 1-irreducible function {v'v')i--„- This case requires 
special attention, sec Rcfs. [54], and from now on we always assume d > 2. Then the inclusion of the counterterms 
is reproduced by the multiplicative renormalization of the action functional (2.1), (2.2) with only two independent 
renormalization constants Zi^2'- 

Sr^ = SyR{v\ v) + 0'Dg0'/2 + 0' [-Vt H- uvZ2d''] 0, (3.2) 

and 

SvR{y', v) = v'DW/2 + v' [-Vt + z/Zi92] V. (3.3) 
In the one-loop approximation the renormalization constants have the forms 

where Sd = Sd/{2T:Y and Sd = 2'k'^I^ /T{d/2) is the surface area of the unit sphere in d-dimensional space. Here and 
below, we use the minimal subtraction (MS) scheme, in which all renormalization constants have the forms "1 4- only 
poles in £." Since the velocity field is not affected by the fields and 0' , the constant Zi is independent of u; the 
one-loop expression (3.4) was presented in [50] and the two-loop calculation was performed much later in Refs. [46,55]. 
The constant Z2 is determined by the 1-irreducible function {0'0)i^ir, which does not involve the correlation function 
(1.3); see item (iii) above. Therefore, Z2 in our model coincides exactly with the corresponding renormalization 
constant for the case of a passive scalar without the random noise, derived in the one-loop approximation in Ref. [56] 
(see also Refs. [57]).^ 



^In the books [16,45], there is a misprint in the expression for Z2 on pages 709 and 115, rcspeetively. It is also interesting 
to note that the one- loop expression for this constant coincides with its analog for a passive magnetic field ("kinetic regime") 
derived in Refs. [58]. 
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Renormalization (3.2), (3.3) can be reproduced by the following multiplicative renormalization of the parameters 
go, uq, and uq: 

go = giJ^^Zg, uo = uZ^, uo = uZu, Zi, = Zi, Z„ = Z2Z^'^, Zg = Z^^. (3.5) 

Here g, u, and v (without a subscript) arc the renormalized analogs of the corresponding bare parameters (with 
the subscript 0) and /x is the reference mass (additional arbitrary parameter of the renormalized theory). The last 
relation in (3.5) is the consequence of the absence of renormalization of the term with Dy in (2.2). The amplitude 
Do in the term with Dy should be expressed in renormalized parameters using the relations Dq = goi/g = gfi'^'^v^. No 
renormalization of the fields $ and "masses" rn = 1/^, 1/L is needed. 

Lot W{e{)) be some correlation function in the original model (2.1) and WR{e, /i) its analog in the renormalized theory 
with action (3.2). Here Cq is the complete set of bare parameters, and e is the set of their renormalized counterparts. 
The relation 5(i>, eo) = Sr{^, e, ji) for the action functional yields H^(eo) = Wji{e, jj) for any correlation function of 
the fields $; the only difference is in the choice of variables and in the form of perturbation theory (in g instead of 
go). We use to denote the differential operation jid^ for fixed eo and operate on both sides of this equation with 
it. This gives the basic RG equation: 

VRGWR{e,n) = Q, VR,G=V^+Pgdg + Mu-l.V,. (3.6) 

Here Vrg is the operation expressed in the renormalized variables, = xd^ for any variable x, and the RG 
functions (the /? functions and the anomalous dimensions 7) are defined as 

7f = In Zf for all F, (3.7) 



fig = 'Df^g = g [-2£ + 37,.] , j3u = V^u = -uju, (3.8) 
where the relations (3.5) have been used. 



IV. FIXED POINT, INFRARED SCALING, AND CRITICAL DIMENSIONS 

It is well known that possible scaling regimes of a renormalizable model are associated with the IR stable fixed 
points of the corresponding RG equations. In our model, the coordinates g^, u* of the fixed points are found from the 
equations 

f3gig*) = f3uig*,u^) = (4.1) 

with the (3 functions from (3.8), while the type of a point is determined by the 2x2 matrix consisting of the elements 
^ = {'-^g0gj duPg, dgPu,du(3u} Calculated at the point g*,u». For IR stable fixed points the matrix Q is positive, i.e., 
the real parts of all its eigenvalues are positive. In our model vanishes identically, and the eigenvalues of the 
matrix CI are simply given by its diagonal elements. 

The analysis of the explicit one-loop expressions shows that, for e > 0, the RG equations for our model possesses 
the only IR attractive fixed point in the physical range of parameters (5*, m* must be positive). The coordinates are 
calculated as series in e, 

g.=gi'h + gPe''+0{e^), = + u^e + ^(e^), (4.2) 

with the one-loop approximation 

(with Sd from (3.4)) for arbitrary d > 2. The two-loop result for was presented in [46,55]: gi^^ ~ — l.Ol.gl^'' for 

d = 3. We have also calculated the two-loop correction for u»: ui^^ « — 0.035wl°'' for d = 3. We shall not expound on 
the derivation of this result because, as we shall see below, the two-loop corrections to the coordinates 5* , are not 
needed for the two-loop calculation of the anomalous exponents. 

Existence of the IR stable fixed point implies that the correlation functions of our model in the IR range exhibit 
scaling behavior with definite critical dimensions Ap of all the fields and parameters F. Let F be some multiplicatively 
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renormalized quantity (a field, parameter or composite operator), that is, F = ZpF^ with certain renormahzation 
constant Zp- Then its critical dimension is given by the expression (see e.g. [16,44,45]) 

A[F] = Af = 4 + A„d^+7^, A„ = 2-7;:, (4.4) 

where d^p and are the corresponding canonical dimensions, 7^. is the value of the anomalous dimension 7f(.9) = 

InZp at the fixed point in question, and A^^ is the critical dimension of frequency. Owing to the exact relation 
between 'j^, and (3g in (3.8), its value at the critical point is found exactly: A^j = 2 — 2e/3 (without corrections of order 

and higher) . As a consequence, the critical dimensions of some basis parameters, fields and composite operators 
are also found exactly: 

A[t;"] = nA„ = n(l - 2e/3), A[r] = nAg = n(-l + e/3), 

A[v'] = d-A,, A[e']=d-Ag, A„ = l. (4.5) 

To avoid possible misunderstanding, it should be noted that simple linear relations for the dimensions of composite 
fields and 6'" follow from the fact that these operators are not renormalized {Zp = 1). For the powers of the 
velocity, this is a consequence of the Galilean symmetry (sc!e [43 45]), while for 9" will be discussed below. 

Let G(r) = {Fi{x)F2{x')) be, for definiteness, some equal-time two-point quantity, for example, the pair correlation 
function of the primary fields $ or some multiplicatively renormalizable composite operators. The existence of a 
nontrivial IR stable fixed point implies that in the IR asymptotic region Ar » 1 and any fixed mr the function G{r) 
takes on the form 

G{r) ~ Uq"^ A'^'^ (Ar)-^° $(rnr). (4.6) 

Here the UV momentum scale A is defined by the relations go = Dq/uq = A^*^, and ^ is certain scaling function whose 
explicit form is not determined by the RG equation itself. The canonical dimensions dg, do and the critical dimension 
Ag ot the function G{r) are equal to the sums of the corresponding dimensions of the quantities Fi^2- 

V. RENORMALIZATION OF RELEVANT COMPOSITE OPERATORS 

In the following, an important role will be played by the composite operators built of the field 0{x) and its spatial 
derivatives. 

We recall that the term "local composite operator" refers to any monomial or polynomial built of the fields $ and 
their derivatives at a single spacetime point x = {t,x}, for example ^" or 9'{vd)6. 

Coincidence of the field arguments in correlation functions containing an operator F gives rise to additional UV 
divergences, removed by a special renormahzation procedure. Owing to the renormahzation, the critical dimension 
Ap associated with certain operator F is not in general equal to the simple sum of critical dimensions of the fields 
and derivatives entering into F. As a rule, composite operators "mix" in renormahzation, that is, an UV finite 
renormalized operator is a linear combination of unrenormalized operators, and vice versa. 

In general, counterterms to a given operator F are determined by all possible 1-irreducible Green functions with 
one operator F and arbitrary number of primary fields, F = {F{x)^{xi) . . . $(a;2))i-ir- The total canonical dimension 
(formal index of divergence) for such quantities is given by 

dr = dp -^Nis,d^, (5.1) 

with the summation over all types of fields entering into the function and the canonical dimensions from Table I. For 
superficially divergent diagrams, dr is a non-negative integer. 

Consider the simplest operators of the form 6"{x) with the canonical dimension dp = —n, entering into the structure 
functions (1.1). From Table I and Eq. (5.1) we obtain dp — —n + Ng — {d — l)Ny> — — {d + l)Ng/, and from the 
analysis of the diagrams it follows that the total number of the fields 6 entering into the function F can never exceed 
the number of the fields 9 in the operator ^" itself: Ng < n {& consequence of the linearity of the original stochastic 
equations in 9). Therefore, the divergence can only exist in the functions with Ny = Nyf = Ng> = 0, and arbitrary 
value oi n — Ng, for which the formal index vanishes, dr = 0. However, at least one of Ng external "tails" of the 
field 9 is attached to a vertex 9'{vd)9 (it is impossible to construct nontrivial, superficially divergent diagram of the 
desired type with all the external tails attached to the vertex F), at least one derivative d appears as an extra factor 
in the diagram, and, consequently, the real index of divergence is necessarily negative. 
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This means that the operator 6" requires no counterterms at all, that is, it is in fact UV finite: 9"^ = Z [O'"-]^' with 
Z = I. It then follows that the critical dimension of 9^{x) is simply given by the expression (4.4) with no correction 
from 7^;, and therefore reduces to the sum of the critical dimensions of the factors: A[^"] = nAg = n(— 1 + e/3), as 
already stated in Eq. (4.5). 

Now let us turn to the scalar operators 

F,, = {dmer (5.2) 

with dp = 0, dp = —n. As we shall see below, it is their critical dimensions that determine the anomalous exponents 
for the structure functions (1.1) and other quantities. In this case, from Table I and Eq. (5.1) we find dp = 

Ng — Nv — (rf — l)Nvi — (d + l)Ng/, with the necessary condition Ng < 2n following from the structure of the 
diagrams. It is also obvious from the analysis of the diagrams that the counterterms to these operators can involve 
the fields 9, 9' only in the form of derivatives, 89, 89', so that the real index of divergence rfp has the form dp = 
d-p — Ng — Ngi — Ny' = —Ny — (d + 2)Ng' — dNy' . It then follows that superficial divergences exist only in the 
correlation functions with = Ny> = Ngi = and any Ng < 2n, and the corresponding operator counterterms are 
reduced to the form Fk with k < n. Therefore, the operators F„ can mix only with each other in renormalization, the 
corresponding infinite renormalization matrix Zp = {Znk} is in fact triangular, Z„k = for k > n, and the critical 
dimensions associated with the operators F„ are determined by the diagonal elements Z„ = Z„„ from the equation 
(4.4) with the anomalous dimension 7„ = C^lnZn. 
Finally, consider irreducible Z-th rank operators of the form 

Fni=8i,9---8i,e{di98i9Y + ■■■, n = l + 2p (5.3) 

with dp = 0, dp = —n/2 (note that i^„o = Fp). Here the dots stand for the appropriate subtractions involving 
the Kronecker 6 symbols, which ensure that the resulting expressions are traceless with respect to any given pair of 
indices, for example, di9dj9 — Sijdk9dk9/d. Of course, the numbers n and I have the same parity, that is, they can 
only be simultaneously even or odd. Like for the operators (5.2), one can show that the operators (5.3) mix only 
with each other in renormalization, the corresponding renormalization matrix is triangular, the critical dimensions 
are determined by its diagonal elements Zni, and the anomalous dimensions are jni = 2'/iln.^„i. 

One important remark is relevant here. The matrix elements Z„ for the operators Fn and Zni for Fni are determined 
by the 1-irreducible correlation functions {F9 ■ ■ ■ 0)i_ir, in which the number of the fields 9 equals to their number in 
the operator F. It is easily seen that the corresponding Feynman diagrams do not involve the bare propagator (99) o 
from (2.6), and, hence, the correlation function of the scalar random noise (1.3). As a result, the critical dimensions 
of the operators (5.2) and (5.3) are completely independent of the form of the scalar forcing in Eq. (1.2). We also 
note that for the same reason the operators (5.3) with equal n and different I do not mix in renormalization: this 
is forbidden by SO{d) symmetry, which is present in the relevant diagrams even if the correlation function (1.3) is 
anisotropic. 

In contrast to (4.5), the critical dimensions A„ and A„; of the operators (5.2) and (5.3) are nontrivial; they are 
calculated as series in e: 



E-'^i'^ (5.4) 



■ 

(of course, Ap = A„o with n = 2p). An important exception is provided by the dimension of the operator Fi = 
{di9di9), the local dissipation rate of the scalar fluctuations, which can be found exactly; see below. The calculation 
to the order (two-loop approximation) will be presented in Sec. VII in detail, and here we give only the first-order 
result: 

(1) ^ n(n-2) {d+l)l{d + l-2) 

3(d + 2) 3(d-l)(d + 2) ■ ^ ' ' 

Expression (5.5) was already presented in [29]. It differs only by an overall factor from its analog for Kraichnan's 
model [9,29] or Gaussian model with finite correlation time [29].^ 



^Morc precisely, the first-order result for A„; in Kraichnan's model is obtained from (5.4), (5.5) after the substitution e 3^/2, 
where g is the exponent in the velocity-velocity correlation function (vv) oc 5{t — t')/k'^~^^. Thus for the "physical" values (e = 2 
and ^ = 4/3) they coincide. 
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The result Ai = A20 = in (5.5) is in fact valid to all orders in e. This can be demonstrated using the Schwinger 
equation of the form 

j [e{x) exp Sr{^) + /5e'{x) = 0, (5.6) 

where Sb. is the renormalized action (3.2) and the notation introduced in (2.3) is used. (We recall that in the general 
sense of the term, Schwinger equations are any relations stating that any functional integral of a total variational 
derivative vanishes; see e.g. Ref. [59].) Equation (5.6) can be rewritten in the form 

{{6' Dee - Wt[e^/2] + uZuZ,d^[0y2] - uZuZ,Fi))a = -Ae'SWR{A)/5Ae, (5.7) 

where Dg is the correlation function (1.3), ((. . .))a denotes the averaging with the weight exp[S'/j($) + A<1>], Wn is 
determined by the Eq. (2.3) with the replacement S Sn, and the argument x common to all the quantities is 
omitted. 

The quantity {{F))a is the generating fimctional of the correlation fmictions with one operator F and any number 
of the fields $, therefore the UV finiteness of the operator F is equivalent to the finiteness of the functional {{F))a- 
The quantity on the right-hand side of (5.7) is UV finite (a derivative of the renormalized functional with respect to 
finite argument), and so is the operator in its left-hand side. Our operator Fi does not admix in renormalization to 
the operator O'DgO {Fi contains too many fields 9) and to the operators Vt[6'^/2] and 9^[6'^/2] (they have the form of 
total derivatives, and Fi does not reduce to such form). On the other hand, the operator 9' DgO does not admix to Fi 
(it is nonlocal, and Fx is local), while the derivatives V^[^^^/2] and 9^[^^/2] do not admix to Fx owing to the fact that 
each field 6 enters in the counterterms of the operators F„ only in the form of derivative d9 (see above). Therefore, 
all three types of operators entering into the left-hand side of Eq. (5.7) are independent, and they must be UV finite 
separately. 

Since the operator vZuZ^Fi is UV finite, it coincides with its finite part, i.e., uZ^Z^Fi = i^F^-, which along with 
the relation Fi = ZiFf gives Z^^ = ZuZy and therefore 71 = — 7„ — 7,^. We recall that at the fixed point 7* = 2c/3 
and 7^ = (the latter equality follows from the relation = — U7„ and > 0), so that 7^ = — 2£/3. From the 
relations (4.4) and Table I one obtains A„ = 2n£r/3 -I- 7*. Combining these expressions gives the desired exact result 
Ai = 0. It will be used later to prove the absence of anomaly for the second-order structure function. What is more, 
it can be used to essentially simplify the calculation of the other dimensions A„i to order e^; see Sec. VII. 



VI. OPERATOR PRODUCT EXPANSION AND COMPOSITE OPERATORS 

The representation (4.6) describes the behavior of the correlation function G(r) for Ar 3> 1 and any fixed value 
of mr. In particular, for the structure functions (1.1) using the data from Table I, Eq. (4.5) and the definitions 
qq = Dq/uq = A^^ one obtains 

5„(r) =Do-"/V"^"en(wO- (6.1) 

For our model, odd structure functions vanish, but they become nontrivial if, for example, the random force in Eq. 
(1.2) is replaced by an imposed constant gradient. The inertial range corresponds to the additional condition that 
mr <SC 1. The form of the function ^(mr) is not determined by the RG equations themselves; in the theory of critical 
phenomena, its behavior for mr — > is studied using the well-known Wilson operator product expansion; see e.g. 
Refs. [15,16]. This technique is also applicable in the theory of turbulence [16,43-45]. 

According to the OPE, the equal-time random quantity in the left-hand side of Eq. (1.1) at x = (xi-|-X2)/2 = const 
and r = xi — X2 ^ can be represented in the form 

[0{t, xi) - 0{t, X2)]" = J2 CF{r)F{t, x), (6.2) 

F 

where the functions Cp are the Wilson coefficients regular in and F are various composite operators (more 
precisely, sec below). 

In general, the operators F entering the expansion like (6.2) are all possible renormalized local operators, allowed 
by the symmetry of the model and the quantity in the left-hand side. In practice, they can be found as the monomials 

which appear in the corresponding Taylor expansions, and all possible operators that admix to them in renormalization. 
If these operators have additional vector indices, they are contracted with the corresponding indices of the coefficients 
Cf. 
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With no loss of generality it can be assumed that the expansion in Eq. (6.2) is made in the operators with definite 
critical dimensions Ap- The structure function (1.1) in rcnormalized variables is obtained by averaging Eq. (6.2) with 
the weight cxpSn with Su from Eq. (3.2), then the quantities (F) appear on the right hand side. Their asymptotic 
behavior for m ^ is found from the corresponding RG equations and has the form {F) oc w'^*'. 

Combining the operator product expansion (6.2) with the asymptotic representation (6.1) we therefore find the 
following expression for the scaling functions ^ni'mi') in the region mr <C 1: 

Umr)=Y,^nF{mr)'^^, (6.3) 

F 

with the coefficients A„F{mr) regular in (mr)'^. 

Obviously, the leading term of the asymptotic behavior of the function (6.3) for (mr) <C 1 is determined by the 
operator with the minimal dimension Ap. The following additional considerations should be taken into account. 

(i) With no loss of generality, it can be assumed that the expansion (6.2) is made in irreducible traceless tensor 
composite operators. In the isotropic case, the mean values of all nonscalar irreducible operators vanish, and their 
dimensions do not appear in the right-hand side of Eq. (6.3). 

(ii) Owing to Galilean invariance of the model and the structure functions (1.1), only invariant operators appear in 
the expansion (6.2). 

(iii) The action hmctional (2.1) and the hmctions (1.1) arc invariant with respect to the shift 9^9 + const, and 
the operators on the right-hand side of (6.2) must also be invariant. This means that they can involve the fields 9 
only in the form of (covariant) derivatives di9 or Vt9. 

(iv) Using the linearity of the equation (1.2) in the field 9, one can show that for any operator F that appear in 
the expansion like (6.2), the number of the fields 9 cannot exceed their total number on the left-hand side. 

Finally, we recall that Ap = dp + 0{e). Thus we may conclude that, at least for small e, the leading terms in the 
small-(mr) behavior for the even function S2n is given by one of the operators Fk from (5.2) with k < n. Indeed, 
any additional derivative or a field different from leads to an increase of the dimension Ap; the operators Fk with 
k > n are forbidden by the property (iv), while the operators containing more fields than derivatives are forbidden 
by (iii). From the explicit form (5.5) it follows that the dimension monotonically decreases as k grows. We finally 
conclude that the leading term in (6.3) is given by the contribution of the operator F^ from (5.2) and 

S2n = (wr)^" a r"+T" (6.4) 

with the dimension 7* defined above Eq. (5.3). 

This establishes the existence of anomalous scaling for the passive scalar field in our model with the identification 
(2n = —2nA0 + A„ = n -I- 7*; see text below Eq. (1.1). 

If the large-scale anisotropy is introduced to the problem by the correlation function of the scalar noise (1.3), tensor 
operators acquire nonzero mean values and their dimensions also appear on the right-hand side of the expansion (6.3). 
An l-th rank irreducible operator gives rise to a term in ^ (mr) proportional to the spherical harmonics Yim for d = 3 
or their analogs for general d. In the special case of uniaxial anisotropy, when the function C in Eq. (1.3) depends on 
a fixed unit vector n in addition to r, they reduce to terms proportional to the Gegenbauer polynomials Pi (Legendre 
polynomials for d = S). 

Repeating the above analysis we conclude that the leading term in the l-th anisotropic sector of the scaling function 
^(mr) for mr <C 1 is determined by the Z-th rank tensor operator Fni (5.3) with the dimension A„( from (5.5). In 
particular, for the case of uniaxial anisotropy 

5„ oc ••• + Pi(cosi?)r"+Tn' +•••, (6.5) 

where 7*; is defined in Sec. V below Eq. (5.3), ^ is the angle between the directions n and r, and the dots stand for 
the contributions of the other anisotropic sectors. It remains to note that the odd functions <S'2n-i-i are nontrivial if, 
for example, the random force in Eq. (1.2) is replaced by an imposed constant gradient, and their leading terms are 
then determined by the vector operators i^2n-i-i,i- 

VII. CALCULATION OF THE CRITICAL DIMENSIONS OF OPERATORS F„i IN THE TWO-LOOP 

APPROXIMATION 

A. General scheme and the relevant diagrams 

From now on, we shall consider composite operators (5.3) in the model without the scalar noise in Eq. (1.2), that 
is, with Dg = in the action functional (2.1). The stirring force in Eq. (1.4), that is, the term with Dy in the action 
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functional (3.2), should be retained. Then the operators (5.3) are renormalized multiplicativcly, Fni = ZniF^^: see 
Sec. V. The renormalization constants Z„( = Zni{g,u,d,s) are determined by the requirement that the 1-irreducible 
correlation function 

{F^i{x)e{xi) . . .^(x„))^_.^. = Z-i'{FMe{xi) . ..e{xn)),_,^ ^ Z-lVni{x;x^, ...,xn) (7.1) 

be UV finite in renormalized theory (3.2), (3.3), i.e., have no poles in e when expressed in renormalized variables 
(3.5). This is equivalent to the UV finiteness of the product Z~^Tni{x; 6), in which 

r„((a;;6') = ^ j dxi . . . j dXnTni{x;xi, . . . ,Xn) 9{xi) . . .9{xn) (7.2) 

is a functional of the field 9{x). In the zeroth approximation, the functional (7.2) coincides with the operator Fni{x), 
and in higher orders the kernel r„; (x; Xi , . . . , a;„) is given by the sum of diagrams shown in Fig. 1 up to the two- loop 
order with their symmetry coefficients (coefficients of the diagrams Nos 2, 4 and 6 are equal to unity, while coefficients 
of the diagrams Nos 8 and 9 arc nontrivial, but they arc not shown because we will not need them; sec below). 
The dashed lines denote the propagators (2.4), while the solid lines denote the propagators (2.6); the slashed ends 
correspond to the auxiliary fields v' and 6' , the ends without a slash correspond to the fields v and 6. Since we are 
working with the renormalized theory, the replacements — s- i^, — > uv should be made in the denominators of 
Eqs. (2.4), (2.6), and the amplitude in (ot)o should be expressed in renormalized variables; Dq = gy^'^v^ (see text 
below Eq. (3.5)). The diagrams Nos 1-7 involve only the vertex (2.7) while the diagrams Nos 8 and 9 also involve 
the vertex (2.5). One dashed line attached to any of the vertices (2.5) must be slashed; there are two variants for the 
diagram No 8 and three variants for No 9. We do not show these variants explicitly (and do not show the slashes and 
symmetry coefficients for these diagrams) , because, as we shall see below, the total contribution of the diagrams Nos 
2, 6, 8, and 9 can be found without the practical calculation of their individual contributions. 

The thick dots in the diagrams correspond to the vertices of the composite operator F^i from (5.3). According to 
the general rules of the universal diagrammatic technique (sec e.g. Rcf. [59]), for any composite operator -^(a;) built 
of the fields 0, the vertex with > attached lines corresponds to the vertex factor 

Vk{x; xu . . . ,Xk) = 6''F{x)/59{xi) . . .Se{xk). (7.3) 

The arguments xi . . . Xfc of the quantity (7.3) are contracted with the arguments of the upper ends of the lines {99')o 
attached to the vertex. For our operators (5.3), built solely of the gradients Wi{x) = d9{x) / dxi at a single spacetime 
point X, the factors (7.3) contain the product di^5{x — Xi) . . . di^d{x — Xk), and the integrations over xi . . .Xk are 
easily performed: the derivatives move onto the upper ends of the corresponding lines {99')o attached to the vertex, 
and their arguments xi . . .Xk are substituted with x. After the derivatives have been moved inside the diagram, the 
remaining vertex factor for the operator F{x) can be understood as a usual derivative: 

Vi,...i, (x) = d>'F{x)/dwi, (x) . . . dwi, (x). (7.4) 

The analysis of the diagrams shows that for any argument Xg in the quantity (7.2), the corresponding spatial 

derivative is isolated as an external factor from each diagram. Using the integration by parts, these derivatives can be 
moved onto the corresponding fields 9{xs), so that the quantity (7.2) can be represented as a functional of the vector 
field Wi = di9: 

^nl{x;9) = ^^ j dxi...j dXnf''^i--''''{x\Xi,...,Xn)Wi^(xi)...Wi^{Xn). (7.5) 

The diagrams that determine the kernel V in (7.5) contain only logarithmic UV divergencies. Therefore, in order to 
find the constant Z~i it is sufficient to calculate the functional F with some special choice of its functional argument 
Wi, namely, one can replace it by its value at the fixed point x, the argument of the operator F^i in Eqs. (7.1). 
Now the product Wi^ (x) . . . Wi^ (x) can be taken outside the integrals over xi, . . . , x„ in (7.5), so that the functional 
F„;(.t; 6) turns to a local composite operator. The integration of the remaining function F„; over .xi, . . . , x„ gives a 
quantity independent of any coordinate variables, and its vector indices can only be those of Kronecker delta symbols. 
Their contraction with the indices of the product Wi^ (x) . . . Wi^ (x) gives rise to the original operator Fni (x) with 
some scalar coefficient F. The integration over xi, . . . , x„ means that in the Fourier representation, the corresponding 
correlation function is calculated with all its momenta set equal to zero, which is always implied in what follows. The 
IR regularization is provided by the parameter m in the function (1.5). In a compact notation one can write: 

Zr,i = ZF, Fni{x) = F, Tni{x;9) = T = FT. (7.6) 
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In the MS scheme all renormalization constants have the form 

■oc- 

Zp = \^Y.^fe-\ (7.7) 
fe=i 

where the coefficients z'^p' = Zp\g, u, d) are independent of s. Then for the corresponding anomalous dimension jp 
from the definitions (3.7) and relations (3.8) for the /3 functions one obtains -jp = V^ln Zp = [(3gdg + /3„9„]lnZir = 

—2T)gZp'' + terms containing poles in e (wc recall that Vg = gdg). However, the pole parts must cancel each other 
owing to the UV finiteness of the anomalous dimension ^p. We therefore arrive at the expression 

= -2VgZ^^\ (7.8) 

that is, in order to find the dimension ^p it is sufficient to find the first-order residue Zp'' in the expansion (7.7). If 

desired, the higher-order residues Zp with k >2 can be calculated to check the aforementioned cancellation (and 
thus the correctness of the calculations). 

Now we turn to the practical calculation of the diagrams needed to determine the coefficients Z^^i in the constants 
Zni from (7.1) to order g'^ (two-loop approximation). 



B. Scalarization of the diagrams and contractions of basic tensor structures 

The contribution of a specific diagram into the functional F in (7.5) for any composite operator F, built of the 
gradients Wi = di6, is represented in the form 

r = Vi2... nt: waWb (7.9) 

where V12... is the vertex factor (7.4), /fl.;; is the "internal block" of the diagram with free indices, the product 
WaUJh ■ ■ ■ corresponds to external lines. The muncrical indices 1,2, .. . will always be imdcrstood as ii,*2, • • their 
number in Eq. (7.9) equals to the number of the letter indices a, 6, . . . and is determined by the number of "rays," 
that is, the number k of lines that attach to the vertex of the operator. These lines are given by products of the 
propagators {Od')o from Eq. (2.6) and are connected by the lines {vv)o and {vv')q from Eq. (2.4); see examples in Fig. 
1. For the two-loop calculation, it is sufficient to consider diagrams with k = 2 and 3, because the diagram No 5 with 
k = 4 factorizes into products of the blocks with k = 2 and, therefore, gives no contribution to the first-order residue 
Zp\ [This is true only if the IR regularization in the correlator (1.5) is provided by the sharp cutoff at fc = m = 1/^, 
so that the one-loop integral is a pure pole in e; see Eq. (7.21) below. If the regularization is provided, e.g., by 
the substitution fc^ ^ fc^ + in (1.5), the one-loop integral contains an 0(e°) term, the diagram No 5 contains a 
first-order pole in e and contributes to ZpK However, the total value of Zp^ in the MS scheme is independent of the 
choice of IR regularization.] 

Since the vertex factor (7.4) and the product WaWb ■ ■ ■ are symmetrical with respect to any permutations of their 
indices, the quantity " in Eq. (7.9) is automatically symmetrized with respect to any permutations of the letter 
indices a,b,. . . and the numerical indices 1,2,.... In what follows, such symmetrization will be denoted by the symbol 
<S. For any fixed number of rays k, the quantity <S / is represented as a linear combination 

SI = J2BiSi (7.10) 

i 

of certain basis tensor structures Si = {Si)i2.'.'. with certain numerical coefficients Bi. There are two such structures 
for the fc-ray diagrams with k = 2 and 3; they have the forms: 

k = 2: Si^S[6ia62bl S2 = S [5i2Sab]; (7.11a) 

k = S: Si = S [dlad2bS3c], S2 = S [Sl2SabS3c]. (7.11b) 

The quantities which will be directly calculated from the diagrams are not the coefficients themselves, but the 
following scalar quantities related to them: 

Ai = tr [{Sirit-SI^t:] = tr [Si -SI], (7.12) 
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where the symbol tr denotes the contraction with respect to aU repeated indices, which will not be shown explicitly. It 
is therefore necessary to express the coefficients Bi in (7.10) in terms of the quantities (7.12). We omit the derivation, 
which is identical to the case of Kraichnan model (see Ref. [19] for detail) and give only the answer: 

k = 2: Bi = 2a2[dAi - A2], B2 = Q!2[-2^i + (d + 1)^2] with 02 = [{d - l)d{d + 2)]-'^; (7.13a) 
k = 3: Bi=6a3[{d + 2)Ai-SA2], B2 = 9a3[-2Ai + (rf + 1)^2] with 03 = [(rf - l)d(rf + 2)(rf + 4)]"^ (7.13b) 

The next step is the contraction of the quantities /fl.;; ^1- (^-9) with external factors: the vertex factor V12... of 
the composite operator and the product WaWb . . .. Again, the derivation is identical to the case of Kraichnan model 
(see [19] for details) and we only give the result: 

r = FT, T = ^kiBi, (7.14) 

where we use the notation (7.6) and the coefficients ki have the forms: 

k = 2 : ki=n{n—l), fc2 = A„;, where Xni — {n — l){d + n + I — 2); (7.15a) 
fc = 3: fci = n(n- l)(n- 2), fca = (n - 2)A„,. (7.15b) 

Finally, combining Eqs. (7.13) and (7.15) we express the function F in the scalar quantities A^: 

T = Y,PiAi, (7.16) 

i 

where 

k = 2: pi=2a2[n{n-2)id-l) + Xi], p2 = a2[n{n + d){d - 1) - {d + l)Xi] ; (7.17a) 

k = 3: pi = 6a3(n- 2) [n(n- 4)(d- 1) + 3A;] , p2 = da^in - 2) [n{n + d){d - 1) - {d + l)Xi] (7.17b) 

with Ai^2 from (7.12), a2,3 from (7.13), and A; = l{l + d — 2). Note that the coefficients pi in (7.17b) vanish for n = 2 
(in general, the diagrams with k rays give no contribution to the functions F for the operators F„; with n < k). Also 
note that the coefficient pi in (7.17a) vanishes for n = 2, Z = 0; this fact will be important in what follows. 



C. General relations for the anomalous dimensions 



Let us denote by J? = Zp^ the first-order coefficient in the expansion (7.7) for the renormalization constant Zp. In 

the perturbation theory, it is calculated as the series R = 9^^s in powers of the renormalized coupling constant 

g, with coefficients Rs depending on u and d. In their turn, they can be written as the sums Rs — X]fe=2 ^s'^ where 
Ri''^ is the total contribution of the diagrams with k rays. As illustrated by Fig. 1, the number of terms in the latter 
sum is always finite: for s = 1, there is the only contribution with = 2, for s = 2, there are contributions with k = 2 
and 3, and so on. Similar decompositions can be written for the coefficients in front of the l/e contributions to the 
scalar quantities A, in (7.12). The corresponding coefficients will be denoted by a^^\ where z = 1, 2 is the number of 
structure, s is the order in g, and k is the number of rays. Prom the definition of the constants Zp and expressions 
(7.6) and (7.9) it follows that, at the same time, a[f^ is the total contribution of the ?'-th structure to the quantity 



Ri^'K that is, = J2iP\ 



■ '^'a-^''. Here p'^''^ are the coefficients (7.17) in which the numbers of rays are explicitly 
shown. 

We shall calculate the first Ri and the second R2 terms (two-loop approximation), which involves diagrams with 2 
and 3 rays. Using the above definitions, expressions (7.8) and (7.17a), and neglecting the terms of order and higher, 
we can write the following representation to the two-ray contribution 7]^ ' to the anomalous dimension jp = 7„; at 
the fixed point (4.1): 



(2) 

Tp' 



'{pi 



+ P2 



g^af^ + 2glaf^ 



} 



with the coefficients pi^2 from (7.17a). Similarly, for the three-ray contribution one can write 

iP=-^gl{p,af^+P2a^^} 



(7.18a) 



(7.18b) 
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with pi^2 from (7.17b). Since diagrams with three rays appear only in order g"^, the contribution of order g in the 

latter expression is absent. We recall that the quantities a-*^ depend on u and d. In (7.18), the substitution u = u* 
is implied. 

The quantities (7.18) should be calculated up to the order s^. The contribution is of order g'^, so that in (7.18b) 
it is sufficient to take the coordinates g,, u» of the fixed point in the lowest-order approximation: g» = gi^'^e, — mI"^ 
with gi^\ ul"'' from (4.3). [We recall that the upper indices for g^ and denote the orders of the expansion in e; 
see Eq. (4.2).] The contribution 7^-' contains terms of order g and g"^. Therefore, in (7.18a) one should take into 
account the leading correction terms to the coordinates of the fixed point, denoted as gi^^ and u^^ in (4.2). We are 
going to show, however, that the quantity (7.18a) can in fact be calculated without knowledge of the coefficients gi^\ 
ui^\ a'21 and 022^. 

We recall that the dimension 7*, for n = 2, Z = is known exactly: 7^" = 739 = — 2e/3; see the end of Sec. V. This 

(2) (3) 

dimension is completely determined by the two-ray contribution from (7.18a), while jp for n = 2 vanishes. The 
coefficient pi in (7.17a) for n = 2, Z = vanishes, while P2 remains nontrivial: P2 = 2/d. Thus from Eq. (7.18a) we 
immediately find the following exact answer for the quantity in the second square brackets: 

= de/6. (7.19) 

This means that the only contribution that survives in the left-hand side is sgi^^a'^i\^^^m, while the 0(£^) contri- 
butions with g^^ , u^^ and 022^ must cancel each other. In order to verify the validity of our calculations, we checked 
this cancellation for d = 3. All the dependence on n and I in Eqs. (7.18) comes from the coefficients pi^2, so that the 
expression (7.19) determines the contribution in the second square brackets in (7.18a) for all n and I. 

Since for n = 2 and 1 = the coefficient pi vanishes, the exact answer for 7|q gives no information about the 
quantity in the first square brackets. However, from the explicit expression for the only one-loop diagram in Fig. 1 

(2) 

it is not difficult to see that the corresponding structure Ai vanishes identically, so that a\i = in (7.18a). Indeed, 
in the quantity Ai the upper (letter) indices of the tensor " in Eq. (7.9) are contracted with its lower (number) 
indices. In the one-loop diagram this leads to the contraction of the momenta at the vertex (7.3) with the transverse 
projector in the correlator {vv)o from (2.4), which depends on the same momentum: PijCk^ki = 0. 

Therefore, the quantity in the first square brackets appears in fact 0{g'^) and, like for jp\ the coordinates 5* and 
It* should be substituted into it only in the leading-order approximations (4.3). 

It remains to note that for the diagrams Nos 2, 6, 8, and 9 the structures Ai also vanish; the mechanism is the 
same as for the one-loop diagram. Therefore, there is no need to calculate these diagrams at all: their nonvanishing 
contributions A2 are known exactly from (7.19) without practical calculation. 



D. Calculation of the scalar quantities Ai 

We shall not discuss the calculation of the scalar quantities Ai from Eq. (7.12) for the all diagrams from Fig. 1 in 
detail, because this definitely would exceed the readers' patience, and give only examples and general ideas. It has 
much in common with the analogous calculation for Kraichnan's model, discussed in Ref. [19] up to the three-loop 
level in great detail. The present calculation differs from that of [19] in a few respects: 

(i) Diagrams Nos 8 and 9 involve the propagators {v'v)o from (2.4) and the vertex v'vv from (2.5); they are absent 
for the case of a Gaussian velocity field (including, of course, the case of Kraichnan's model). 

(ii) Diagrams Nos 6 and 7 are present for any Gaussian velocity field with finite correlation time. However, for 
Kraichnan's model they effectively contain closed circuits of retarded propagators {0'9)o and therefore vanish. It is 
crucial here that for Kraichnan's model the velocity correlator contains the S function in time. In our case, the velocity 
has finite correlation time and these diagrams give nonvanishing contributions. 

(iii) In Kraichnan's case, the diagram No 2 and, in general, all diagrams with the self-energy insertions are easily 
taken into account: it is sufficient just to drop them, and in the remaining diagrams to replace the bare viscosity i/q 
with its exact analog; see Ref. [29]. This is a consequence of the fact that, in Kraichnan's case, the self-energy (that 
is, nontrivial part of the 1-irreducible function (6''0)i_ir) is exactly given by the simplest one-loop diagram (the other 
contain closed circuits of retarded propagators and vanish). That diagram does not depend on its frequency argument 
and is simply proportional to p^, its squared momentum argument, and thus its contribution only leads to a certain 
redefinition of the viscosity coefficient. 
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In the case at hand, the one- loop self-energy diagram is a nontrivial function of p and the calculation of such diagrams 
also becomes nontrivial. In higher orders, diagrams with multiloop self-energy insertions (absent for Kraichnan's case) 
will also appear. 

(iv) Diagrams Nos 3 and 4 are present and nontrivial both for our model and for Kraichnan's case (as well as for 
the Gaussian model with finite correlation time). Of course, the corresponding analytical expressions are different in 
these two cases due to the difference in the explicit forms of the correlation functions {v'v)o. In particular, all their 
contributions for the zero correlation time are expressed in terms of hypergeometric functions (see Ref. [13]) while for 
our case this is no longer true (see below). 

(v) In Kraichnan's model, the value of was given by the one-loop approximation exactly; in our case, the 
higher-order contributions are nontrivial and should be taken into account. 

As we have seen in Sec. VII C, in the two-loop approximation the total effect of the diagrams Nos 2, 6, 8, and 9 
and of the 0(£^) contribution in 5* can be foimd from the exact identity (7.19) without the practical computation, 
but in the three- loop and higher orders the items (i), (iii) and (v) become nontrivial. 

Consider the calculation of the one- loop diagram. No 1 in Fig. 1. Using the explicit forms of the propagators and 
vertices in (2.6), (2.7) and the definition we obtain (in renormalized variables) 

where the factor 1/2 in front of the integral is the symmetry coefficient (see Fig. 1) and the three factors in the 
integrand come from the vertex of the composite operator, the propagator {vv)o of the velocity field, and the product 
of two propagators {0'6)o, respectively. The second equality in (7.20) is the result of elementary integration over w. 
Obviously, Ai = /"^ = 0, the fact already mentioned in Sec. VII C, while for A2 = I^^ one obtains 

_ ffM'-5rf(d-l) dk _ gSd{ii/mf-{d-l) 

^ Au{u+l) y„ 8u(u+l)£ ^ ■ ' 

with Sd from (3.4); the pole part of this expression is simply obtained by the replacement {p/mY'^ 1. Substituting 
this expression into (7.18a) along with the expressions (4.3) for the fixed points gives the one-loop result (5.5). 

Now let us turn to the calculation of the two- loop diagrams Nos 2-7 in Fig. 1. These depend on two integration 
momenta k and q and two frequencies, which can always be assigned to the two (f w)o lines of a diagram. The integra- 
tions over the frequencies (or, equivalently, the times in the time-momentum representation) are always elementary. 
The result can be interpreted as a sum of terms, each of which corresponding to one possible "temporal version" of the 
diagram; different temporal versions correspond to all possible orderings of the integration times in the diagram. To 
each version corresponds an "energy denominator," given by the product of the factors corresponding to all "temporal 
cross-sections" of the diagram; to each cross-section corresponds the sum of "energies" 5k = vk^ for all intersected 
{vv)q or {v'v)q lines and £k = uvk"^ for {96')q lines. Thus, with some experience, it is possible to write down the 
result of the temporal integration without performing the actual integration. As a typical example, we give the result 
of temporal integration for the quantity /fl corresponding to the diagram No 3: 

rab ^ (W-^)- [ f dkdq {k + q)„Xk + q},q,MjP,j(k)I\2{q) ( I 1 1 . . 

12 4u(l + «)(27r)2d J J {kqY+2^ q2{k + q)2 [P + uq^ + u(k + q)2] \ [k^ + q-^ + u{k + q)2] «g2 J ' ^ ■ > 

now the corresponding scalar quantities Ai are easily obtained. It is convenient to represent the denominators as 
products of simpler factors, and to combine the quantities Ai corresponding to different diagrams; this sometimes 
leads to noticeable simplifications of the integrands. With the only exception (see Eq. (7.32) below), all these 
quantities can be reduced to linear combinations of the following "basis" scalar integrals: 

J J {kqY+2^ k2+q2 + 2(3{k<i)- 8e ^pl/^i- l^-^-^J 



(2^) 



2d 



where the parameter /3 takes different values: /3 = 1, /? = u/(l -|- w) or /3 = ('u/2(l -|- ujY^^, while p= 1,2, and 3, and 
■d denotes the angle between the vectors q and k, so that (kq) = kqcos-d. [We do not discuss much simpler integrals, 
e.g. those that can be factorized into two independent integrals over k and q, and so on.] The integrands in (7.23) 
involve three independent parameters, the moduli k and q and the angle 1?, so that the integrals can be written as 



{2n) 



2d 



dk j dciJ^{k,q,'&) = Sl j dkk'^-^ j dqq'^'^ {J'{k,q,d)) 
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where the brackets denote the angular averaging over the unit sphere in d dimensions normalized such that (1) = 1. 
Let us expand the integrands in (7.23) in /3 or, equivalently, in the scalar product (kq) = kq cosd. In each term of 
the resulting expansion, the integrations over the angles can be computed using the following formulas: 

(cos'" ^) = ... J^^'Ti^'l ^. (7-24) 

with n = 1, 2, The remaining integrals over the moduli have the forms: 

Using the identity 

In{m) = -^Vmlnim), = md/dm, (7.26) 

which follows from the last equality in Eq. (7.25), the integral /^(m) can be represented in the form 

Um) = ^l ^(^j , (7.27) 

that is, the number of integrations is reduced and the pole in e is isolated explicitly. We need only the pole part of 
the integral /„(!), which now is simply obtained by setting £ = in the integrand of (7.27). The resulting integral is 
easily calculated: 

Thus we have represented the pole part of the integrals (7.23) as infinite series with known coefiicients. It is not 
difficult to see that these series can be reduced to the hypergeometric function 



ab a(a+ 1)6(6+1) 



z 



2 



Fia, 6; c; .) = 1 + - . + ,(,V\) " ^ + • • • - (7-29) 



namely, for the quantities ^'2p(/3) defined in (7.23) one obtains 



*'^(^) - r{d/2 - l/2)r(d/2 + ^^p^(inhl,d/2 + l+p,P ), (7.30) 



with Euler's F function. For the first special values of p this gives: 



and so on. The integral 



_J_ r r rfkrfq (kg) fc- sin- ^ _ m-^^gj p 

(27r)2rf y y (fc9)'i+2. (k + q)2 [(k + q)2 + g2 + .^p] - 2e y'-"^' 

(where x = 1/u) does not reduce to the hypergeometric function and can only be expressed in the form of a single 
convergent integral, suitable for numerical calculation, for example 

= 2dill'2)id + 4) I! il + z)il + .zr ^''^^'^ + "'^"^ ^'-''^ 
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or 



J{x,d) 



nd/2) 



az 



)d/2 



^T{{d - l)/2) Jo {x - 1)2 + 4x^2 



2^1/2 



In 



1 + X 



— z{x — 1 + 2z^) arcsinz- 



[2(l+a;)-22]i/2 



arctan 



z[2(l + a;)-z2]i/2 

(l+X-^2) 



(7.34) 



It remains to note that in Eq. (4.3) of Ref. [30] the latter integral is given with a misprint. 



VIII. ANOMALOUS EXPONENTS TO ORDER 

Using the techniques described in the preceding Section, we have performed complete two-loop calculation of the 
critical dimensions A„; of the composite operators (5.3) for arbitrary values of n, I, and d and obtained the following 
expression for the second coefficient in expansion (5.4): 



.(2) 



9(d- l)2(d + 2)2(d + 4) 
+ 9C[n(d + n)((i- 1) - \i{d+l)\ 



|2(d + A)A [n{n - 2){d - 1) + A,] + {n- 2)^68 [n{n - 4)(rf - 1) + 3A;] + 



(8.1) 



with Xi = l{d + I — 2) and 



(d+1) 



^3 



2{d+2){l-x) 2(d-|-4)(l -x)a;(l -|-a;)2 •'V(a; + l)V - d)il - x) 



^J^^^^Jix,d), 



B = 
C = 



(d+1) 



3(l-a;)2(rf + 4) 
1 



xHd + 1) (I 

(rf + 4) ■ 



X+l ^ \2{X+1)) (.T + l)2 

Ax{d+l) 



1 



(X + 1)2 
1 



3x2(d-l) n\ x[2d - I + x{d - 2)] 
9(l-a;)2 1^ 4 ix + l) ^\2{x + l) 



[d+l + 2x{d - 2)] ^ 

H 7 r~7T75 ^2 



(X + 1)2 



{x + iy 



1 



4(rf+l) 



^3 



(a; + l)(rf + 4) •'V2(a; + 1); (a; + l)2(d + 4) ^V(a; + 1)2 



^.2) 



Here J{x, d) is the integral (7.34), x = with from (4.3), and Fk{z) = F{1, 1; d/2+k; z) is the hypergeometric 

function (7.29). The values of Fk entering into (8.2) can be related by the recurrent relation 

{z - l)F2{z) = z{d + 2)Fs{z)/{d + 4) - 1, 

but the resulting expressions look more cumbersome and we have kept both F2 and F3 in the formulas. 

Contributions with A, B, and C in (8.1) come from the structures Ai with k = 2, Ai with A: = 3, and A2 with fc = 3, 

(2) 

respectively. The structure A2 with k = 2 gives no contribution to A^^^' , as discussed in Sec. VII C in connection with 
Eq. (7.19). For the most interesting case d = 3 one obtains: 



A = -0.90239, B = -0.135498, C = 0.19622, J = -0.024976. 
Expression (8.1) simplifies for the most important case of the isotropic sector (even n and I = 0): 



. (2) ^ n{n - 2) 

"° (d-l)(d + 2)2(d + 4) 



2{d + A)A + 6(n - 4)B + 9{d + n)C 



For the simplest nontrivial case n = 4 one obtains 

Afo^ = 8(2^ + 9C)/(d - l)(d + 2)2, 



U) 



(8.4) 



(8.5) 



that is, the quantity B does not enter into the result. For n > 6, all the coefficients (8.2) contribute to the result. 
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IX. DISCUSSION AND CONCLUSION 



We have studied a model of a passive scalar field, governed by the diffusion-advection equation (1.2) and subject to 

a large-scale random forcing (1.3). The advccting velocity field obeys the Galilean-invariant Navicr-Stokes equation 
(1.4) subject to an external random force, white in time and having a power-law spectrum oc fc'*"'*"^^; see Eqs. (1.5) 
and (1.6). 

Using the RG and OPE methods, we have shown that the structure functions of the scalar field display anomalous 
scaling behavior; see Eqs. (6.4) and (6.5). The corresponding anomalous exponents A„ are identified with the critical 
(scaling) dimensions of certain composite fields (operators), namely, powers of the local dissipation rate of scalar 
fluctuations (5.2), which offers the possibility to calculate them within a regular perturbation theory, as series in e; 
see Eq. (5.4). 

The calculation has been accomplished to the second order, (two- loop approximation), including the exponents 
of anisotropic contributions (6.5). The latter are identified with the critical dimensions of tensor composite fields 
built of the scalar gradients (5.3). The first-order expressions (5.5) coincide with the exponents of the well-known 
Kraichnan's rapid-change model (up to a simple normalization), while in the second order they are different. Like for 
the rapid-change model, the second-order structure function is not anomalous. 

Thus we have overcome two important limitations of the previous treatments of the problem: absence of time 
correlations and Gaussianity of the advecting velocity field. It is interesting to note that both the RG-mechanism of 
the anomalous scaling and the results for the exponents are, in many respects, similar to the case of the rapid-change 
model. Let us compare our findings with those for the Gaussian models. 

Universcdity: Independence of the forcing and relevance of the zero-modes picture. As we have seen, the 
critical dimensions of all composite operators (5.2) and (5.3), and therefore the corresponding anomalous exponents 
(including anisotropic sectors), are independent of the forcing, specified by the correlator (1.3). In particular, this 
means that they remain unchanged if the stirring noise in Eq. (1.2) is replaced by an imposed constant gradient, like 
e.g. in Refs. [21,29,33]. The role of the forcing is to maintain the steady state of the system and thus to provide 
nonvanishing amplitudes for the power-law terms with those universal exponents. This behavior is already well known 
for the passive scalar or vector fields, advected by the Gaussian velocity fields with vanishing or finite correlation time. 

In the language of the RG (which is equally applicable to the case of a zero or finite correlation time of the advecting 
field) this is explained as follows: the stirring force or the imposed gradient do not enter into the diagrams that 
determine the renormalization of the operators (5.2) and (5.3), so that their dimensions appear forcing-independent. 
Similar diagrams determine the contributions of those operators into the operator-product expansions (6.2), which 
are nontrivial even for the unforced model. The difference is that for the unforced model, the mean values of the 
operators vanish, so that they give no contribution to the right-hand sides of representations like (6.3). For the 
isotropic correlator (1.3), scalar operators accjuire nonzero mean values and contribute to the right-hand side of (6.3), 
while for the anisotropic correlator or the imposed constant gradient, the mean values of irreducible tensor operators 
also become nonzero and their contributions are "activated" in representations (6.3). 

For the case of a Gaussian advecting field with vanishing correlation time, when the equal-time correlations functions 
satisfy exact closed differential equations, the above picture it is easily understood in the language of the zero-mode 
approach [14]: forcing terms do not affect the corresponding differential operators; thus the anomalous exponents, 
determined by the zero modes (solutions of homogeneous unforced equations) also appear forcing- independent. On 
the contrary, the amplitudes are determined by the matching of the inertial-range zero-mode solutions with the forced 
large-scale solutions, which is only possible in the presence of the forcing terms. 

The exact resemblance in the RG picture of the rapid-change models and the finite-correlated cases suggests that 
for the latter, the concept of zero modes (and thus of statistical conservation laws) is also applicable, although the 
corresponding equations arc not differential and involve infinite diagrammatic series. 

Hierarchy of anisotropic contributions. In the presence of large-scale anisotropy (that is, the anisotropy 
introduced at scales of order L by the forcing in Eq. (1.2)), structure functions of the scalar field can be decomposed 
in irreducible representations of the rf-dimensional rotation group SO(rf). Such a decomposition naturally arises from 
the corresponding OPE, provided it is made in irreducible traceless tensor composite operators; the rank I of a tensor 
operator can be used to label the terms of the SO(d)-expansion and can be viewed as the measure of anisotropy of 
the corresponding term ("sector"). Thus each anisotropic sector is characterized by its own set of scaling exponents, 
the leading term is given by the ?-th rank composite operator with minimal critical dimension. 

Explicit expressions for these dimensions, derived to second order in £, exhibit an hierarchy related to the degree 
of anisotropy: the higher is the rank of the operator (the more anisotropic is the contribution), the larger is the 
corresponding dimension, and thus the less important is its contribution to the inertial-range behavior. This hierarchy 
can be expressed by the relation dlS.ni/dl > 0, which is obvious from the first-order expression (5.5). It holds for all 
values of n and d. This picture is similar to the hierarchy relations derived earlier for the passive scalar and magnetic 
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fields advectcd by the Gaussian velocity ensembles [21,29-32]. 

In particular, this means that the overall leading term is given by the exponent from the isotropic sector, and it 
is therefore the same for the isotropic and anisotropic forcing. It also should be stressed that the independence of 
the scaling behavior in different sectors is a direct consequence of the linearity of our model, independence of the 
exponents on the random force, and the SO(d) symmetry of the unforced model. On the contrary, the hierarchy of 
the exponents follows from the explicit expressions, obtained only by practical calculation. 

According to the Kolmogorov-Obukhov theory [1,2], the anisotropy introduced at large scales by the forcing (bound- 
ary conditions, geometry of an obstacle etc) dies out when the energy is transferred down to smaller scales owing 
to the cascade mechanism (isotropization of the developed turbulence in the inertial-range). The analytical results 
discussed above confirm this classical concept and give a more quantitative picture of the isotropization. 

The hierarchical picture, derived here and in Refs. [21,29-32] for passively advected fields, appears unexpectedly 
general, being compatible with that established recently in the field of NS turbulence, on the basis of numerical 
simulations and natural experiments; see Refs. [60] and references therein. There, the velocity structure functions 
were decomposed in the irreducible representations of the rotation group. It was shown that in each sector of the 
decomposition, scaling behavior can be found with apparently universal exponents. The amplitudes of the various 
contributions are nonuniversal, through the dependence on the position in the flow, the local degree of anisotropy and 
inhomogeneity, and so on. 

It is worth recalling here that the so-called "additive fusion rules," hypothesized for the NS turbulence in a number 
of papers (sec e.g. Rcf. [61]) and characteristic of the models with multifractal behavior (see Ref. [62]), arise naturally 
in the context of the rapid-change models owing to their linearity. The existing results for the Burgers turbulence can 
also be interpreted naturally as a consequence of similar fusion rules, where only finite number of dangerous operators 
contributes to each structure function; see Ref. [63]. This is rather surprising because the equations for the correlation 
fimctions in such cases are neither closed nor isotropic and homogeneous. One can thus speculate that the anomalous 
scaling for the genuine turbulence can also appear, in some sense, a linear phenomenon. Of course, one should not 
insist too much on this bold assumption. 

Universality: Independence of the time scales. An important issue is that of the universality of anomalous 
exponents. As already discussed, the exponents A„; = A„;(e, d) in Eq. (5.4) are independent of the forcing in the 
scalar equation (1.2), and thus independent of all the parameters that can appear in its correlation function (1.3). 

However, the exponents depend on the exponent £ that enters the correlation function of the stirring force (1.5) 
in the NS equation (1.4). They also depend on d, the dimensionality of the x space (note that the basis dimensions 
related to the velocity field are d-independent; see Eqs. (4.5)). 

Earlier, it was argued on phenomenological grounds that the anomalous exponents of the scalar field can depend 
on more characteristics of the advecting field than only the exponents; see e.g. the discussion in Ref. [34]. Indeed, 
analytical derivation of the anomalous exponents of the passive scalar field, advectcd by a Gaussian velocity with 
finite correlation time, has revealed for some asymptotic regimes ("local turnover exponent") their dependence on the 
correlation time of the velocity field (more precisely, the dimensionless ratio of the correlation times of the scalar and 
velocity fields); see Refs. [29,30,35]. 

In our case, the exponents could depend, in principle, on the analogous dimensionless parameter wq = kq/vq from 
(3.1), the (inverse) Prandtl number. After the RG resummation, this parameter is replaced with the corresponding 
invariant variable, which has exactly the meaning of the ratio of the scalar and velocity correlation times (for a detailed 
discussion of this point see Ref. [29]). However, the analysis of the RG equation shows that in the IR asymptotic 
range, this parameter tends to a fixed point, whose coordinate depends on d and e, but not on the initial value uq; 
see Eqs. (4.2), (4.3). As a result, all the dimensions, including A„( from (5.4), appear independent of Uq. In the RG 
language, the nonuniversality (that is, the dependence on the ratio uq or its analog) of the exponents in the Gaussian 
model is a consequence of the infinite degeneracy of the IR stable fixed point; see the discussion in [29]. In the NS 
model, the fixed point is unique (non-degenerate), and the exponents appear universal. 

We stress that, although the coordinates of the fixed point are known only in the two- loop approximation (see the 
discussion below Eqs. (4.2), (4.3)), the statement about the universality is exact, that is, it holds to all orders of the 
e expansion. 

Since the degeneracy of the fixed point in the model studied in Refs. [29,30,35] is an artifact of the Gaussianity of 
the velocity ensemble, wc bcilieve that our result for the non-Gaussian velocity ensemble, described by the Galilean- 
covariant NS equation, suggests that for the real passive advection the anomalous exponents are universal, that is, 
independent of the Prandtl number or the ratio of the scalar and Ivelocity correlation times. This is probably the 
most important qualitative conclusion that can be inferred from our analysis. It is then relevant to discuss the role 
played by the Galilean symmetry of our model in the RG-analysis. 

Sweeping efFects and the Galilean invariance. The results obtained within the RG and OPE approach and 
within the e expansion, are reliable and internally consistent for asymptotically small e. A serious question is that of 
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the validity of the e expansion for finite e's, and the possibihty of the extrapolation of those results to the physical 
value e = 2. 

For the rapid-change model, the e expansion works surprisingly well. It was demonstrated [19] that the knowledge 
of three terms allows one to obtain reasonable predictions for finite e ~ 1; even the plain e expansion captures some 
subtle qualitative features of the anomalous exponents established in the exact solutions of the zero-mode equations 
and in numerical simulations. The quantitative agreement can be achieved with the aid of various improvements, like 
the inverse e expansion or interpolation formulas [19]. 

In the case of the Gaussian model with a finite correlation time, however, there is a natural upper bound for the 
range of validity of the e expansion: for e larger than certain threshold value Sc, the velocity field (and hence all its 
powers) become dangerous: their critical dimensions, known exactly due to the Gaussianity, become negative, and new 
strong IR singularities occur in the diagrams; see the discussion in Ref. [29]. This leads to a qualitative changeover in 
the small-mr behavior of the scalar field, as demonstrated in Refs. [37,38,41] using certain nonperturbative analytical 
methods and numerical experiments. Therefore, the results obtained within the plain e expansion no longer apply. 

Physically, this is a manifestation of the fact that above £c, the so-called sweeping effects (kinematic transfer of the 
small-scale turbulent eddies by the large-scale ones) become important. Thus such threshold value can also be viewed 
as the upper bound of the range of validity of the model itself: the lack of Galilean covariance becomes a serious 
shortcoming of the synthetic Gaussian velocity ensemble when the sweeping effects become important. 

In the model (1.4), (1.5), the dimensions of the powers t;" are known exactly, A[t;"] = nA[w] = n(l — 2e/3); they all 
become negative for e > = 3/2 [43-45]. Some operators, built of the velocity field and its temporal derivatives, also 
become dangerous for e < 2; see [44,45]. Their contributions to the OPE for the correlation functions of the velocity 
and scalar fields become singular; however, they can be summed out explicitly using certain infrared perturbation 
theory. This indeed results in a qualitative changeover in the small-mr behavior of the correlation functions, their 
strong dependence on the IR scale i = 1/m, and superexponential decay in time [43-45], in agreement with the 
phenomenological analysis of Refs. [64]. 

Galilean symmetry of our model guarantees, however, that the invariant quantities, for example, the equal-time 
structure functions (1.1), are not affected by the sweeping. More formally, the contributions of the aforementioned 
dangerous operators do not appear in the OPE for Galilean invariant correlation functions; see Refs. [16,43-46] for 
detailed discussion. This means that in model (1.2) (1.6), the scaling relations obtained for small e, for Galilean 
invariant quantities can be extrapolated beyond the threshold £c, in spite of the fact that the sweeping becomes 
important there. The most recent numerical simulations of the model (1.4)-(1.6) have shown that the scaling relations, 
obtained by the RG analysis for the structure functions, remain valid for e as high as e = 7/4 [51]. 

Extrapolation to the physical value e = 2. Relevance of the model for the real turbulent advection. 
Our calculation of the anomalous exponents implied smallness of the RG expansion parameter e. For small e, a 
serious flaw shared by our model with the Gaussian ones is that the advecting velocity field is non-intermittent, in 
contrast to the real turbulent fluid. However, numerical simulations of Refs. [51,52] suggest that, as e increases, the 
behavior of the model (1.4)-(1.6) undergoes a qualitative changeover and the scaling of the velocity field becomes 
anomalous: the exponents of the structure functions become different from the results of naive extrapolation of the 
small-£ prediction. In the RG language this probably means that certain Galilean invariant operators acquire negative 
critical dimensions for some finite values of e, close to the physical value e = 2. Unfortunately, identification of those 
operators and calculation of their dimensions on the basis of the model (1.4)-(1.6) lies beyond the scope of the present 
RG technique: the effect takes place for finite, and not small, values of e, while the dimensions of the operators 
arc known only in the form of the first terms of the expansions in e (some dimensions arc known exactly, but they 
all remain positive for £ < 2). Detailed discussion of the critical dimensions of Galilean invariant operators can be 
found in Refs. [16,44,45] and the original papers [50,65]. Hopefully, the problem will be solved with the aid of an 
alternative perturbation theory (the expansion in l/d seems very promising, but so far it has been constructed only 
for Kraichnan's model and only to the leading order [10]). 

If the dangerous Galilean invariant operators indeed arise in the model (1.4) (1.6) for some finite values of e, they 
will also contribute to the OPE's for the structure functions (1.1) of the scalar field. Physically, this corresponds to 
the contribution of the velocity to the intermittency of the scalar field, while the contributions of the operators (5.2), 
(5.3) correspond to the intrinsic intermittency of the scalar field itself. Obviously, only the latter contribution can be 
described within the £ expansion. Since the scalar fields appear much more intermittent than the velocity field, one 
can assume that this latter contribution dominates the anomalous behavior of the scalar, or, at least, it is relatively 
more important than the former one. One can therefore hope that the dimensions A„ and A„(, taken at the physical 
value £ = 2, can be identified with the leading anomalous exponents of the structure functions of the real passive 
scalar field. 

Experimental results for the structure functions of a passive scalar field are presented in Refs. [4-6] in terms of 
the exponents (n, with S'n oc r '»". For an even function, in our notation C,2n = —^nAg + A„, where A^ = —1 -|- £/3 
from (4.5) is the critical dimension of the scalar field and A„ are the dimensions of the operators (5.2), with the 
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second-order expression given in Eq. (8.4). For the physical value e = 2 this gives (211 = 2n/3 + A„(£ = 2). 

The results of [4-6] seem to be consistent with the Kolmogorov value (2 = 2/3 for the second-order function, 
in agreement with our exact result Ai = 0. Possible deviation, if any, can be attributed to the anisotropy of the 
experimental setup or/and contribution of velocity's intermittency, neglected in our analysis. 

From the Fig. 3 presented in the most recent study [6] one can infer A2 ^ —0.23 and A3 ~ —0.67, which shows 
clear deviation from the Kolmogorov values (A„ = for all n). For e = 2 and d = 3 the one-loop approximation 
(5.5) gives A2 ^ — 1 and A3 ~ —3; the two-loop correction (8.4) appears numerically very small and does not affect 
this result markedly. Admittedly, it is difficult to speak about a good quantitative agreement with the experimental 
values. 

The situation resembles that encountered for Kraichnan's model. There, the one-loop results for the most realistic 
value ^ = 4/3 are also equal to A2 —1 and A3 ~ —3 (see the remark and footnote below Eq. (5.5)) and essentially 
overestimate the real value of the exponents A2 ~ —0.3 and A3 ~ —0.7, known from the numerical simulations of 
Refs. [42]. The second-order corrections (although different from their analogs in our model) also appear too small 
to improve the agreement. A better agreement is achieved if the third-order correction (which is not small) is taken 
into account, and the straightforward expansion in ^ is augmented by additional considerations about its nature, 
convergence properties, character and location of singularities etc; see Ref. [19]. 

In contrast to Kraichnan's model, there is no reason to believe that in our case the series in e have finite radii 
of convergence. As in most field theoretic models, they can only be asymptotical series. Thus, one should not have 
expected that the straightforward summation of the first two (or even more) terms would give a good result. In models 
of critical behavior accurate theoretical predictions for the exponents imply knowledge of the large-order asymptotic 
behavior of the coefficients of the e expansions, obtained using the instanton calculus, and special resummation 
procedures for the divergent series [15]. In dynamical models the corresponding methods are at their infancy [66-68]: 
to the best of our knowledge, an instanton-typc solution for an action functional of the Martin-Siggia-Rose type 
has been obtained only for a model whose equal-time correlation functions correspond to a system in thermodynamic 
equlibrium [68]. The instanton analysis of Refs. [67], performed in Lagrangian variables (which implies zero correlation 
time of the velocity field, that is, only Kraichnan's case) did not touch upon the problem of the large-order coefficients 
of perturbative series; it has mostly been concentrated on the behavior of the exponents (n at large n. One can hope 
that further development of the instanton techniques for dynamical models, combined with the RG framework will 
give the solution of this important problem. This work is left for the future. 
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TABLE I. Canonical dimensions of the fields and parameters in the model (2.1), (2.2). 
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FIG. 1. Diagrams of the function F from (7.2) in the one-loop and two-loop approximations. 
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